Phase diagram and critical endpoint for strongly-interacting quarks 
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We introduce a method based on the chiral susceptibility, which enables one to draw a phase dia- 
gram in the chemical-potential/temperature plane for strongly-interacting quarks whose interactions 
are described by any reasonable gap equation, even if the diagrammatic content of the quark-gluon 
vertex is unknown. We locate a critical endpoint (CEP) at (fj, E ,T E ) ~ (1.0, 0.9)T C , where T c is 
the critical temperature for chiral symmetry restoration at /j, = 0; and find that a domain of phase 
coexistence opens at the CEP whose area increases as a confinement length-scale grows. 
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A central goal of the worldwide programme in rclativis- 
tic heavy ion collisions is to chart the phase diagram of 
QCD in the plane of nonzero temperature (T) and chem- 
ical potential (/i). This will provide fundamental insight 
into the origin of observable mass and the nature of the 
early universe. Two decades of intense speculation have 
led to an expectation that the phase diagram is complex. 
The existence of a critical endpoint in the (/i, T)-planc 
has been conjectured [1]. In the chiral-limit theory, this 
marks the end of a line of second-order chiral-symmetry- 
restoring (and possibly deconfming) transitions, originat- 
ing on the temperature axis, and the beginning of a line 
of first-order transitions. Such a critical endpoint would 
have observable consequences [2]. Hence it is imperative 
to demonstrate its existence, determine its location and 
demarcate the subsequent domain of phase coexistence. 

Attempts have been made using latticc-QCD. Such 
studies rely on Monte Carlo methods but the absence 
of a probability measure at fi ^ precludes direct com- 
putation. Therefore mathematical devices are necessar- 
ily employed in the search for a CEP. They yield [3-6]: 
H E /T C = 1.0-1.4 and T E /T C w 0.93, and a signal for 
a material phase coexistence region [5]. However, it is 
not yet certain whether the existence of a CEP survives 
in simulations with lattice parameters that more closely 
resemble the physical world [7]. 

Models have also been used to search for a CEP. The 
Nambu-Jona-Lasinio type yield [8, 9]: /^ E /T c w 1.7, 
T E /T C « 0.4; their Polyakov-loop extensions produce 
[10-13]: /i E /T c = 1.5-1.8, T E /T C = 0.3-0.8; and a chi- 
ral quark model gives [14] {fi E ,T E )/T c = (2.0,0.4). On 
the other hand, a Polyakov-loop-augmcntcd chiral quark 
model produces [15] (fi E ,T E )/T c = (0.9,0.8). The for- 
mer, mutually consistent results for the CEP's location 
conflict markedly with those obtained from lattice-QCD: 
H E /T C is significantly larger and T E /T C , much smaller. 
If they are nevertheless correct, then finding the CEP in 
experiment will be difficult because modern colliders are 
restricted to exploration of the small- yU domain. Given 
this observation, it is unsurprising that an analysis of 



flow data from the rclativistic heavy ion collider leads to 
the estimate [16]: fi E /T c > 1.0 and T E /T C < 1.0. 

The Dyson-Schwinger equations (DSEs) provide a non- 
perturbative approach to studying continuum-QCD [17] 
and have been used to prove exact results relating to 
chiral symmetry [18-20]. Simple DSE truncations have 
been applied to the CEP problem. A confining zero- 
width momentum-space interaction, the antithesis of the 
NJL-modcl, produces [21] ^ E /T c = 0, T E /T C = 1; and a 
separable-interaction [22]: fi E /T c = 1.09, T E /T C = 0.78. 
However, neither study described a region of coexisting 
phases. Notwithstanding that, in this chain of remarks 
about the model results, there is a hint that the length- 
scale characterising confinement in the quark-antiquark 
interaction markedly influences the location of the CEP. 

Herein we employ the DSEs to produce a phase di- 
agram for strongly-interacting quarks, to locate a CEP 
and demarcate the coexistence region. The basic tools 
are the chiral susceptibility and the gap equation. In 
QCD its kernel is defined by a contraction of the dressed- 
gluon-propagator and -quark-gluon-vertex. For the for- 
mer we use a form that can interpolate between models 
of the non-confining NJL-type and the confining interac- 
tions used in efficacious DSE studies of hadron observ- 
ables [24] , whilst always providing a super- renormalisable 
interaction. For the latter, we use either the rainbow- 
truncation; i.e., the leading-order term in a symmetry- 
preserving scheme [23], or a dressed- vertex Ansatz. The 
capacity to draw the phase diagram derived from an ar- 
bitrary dressed- vertex is essentially new. 

At T =/= =/= /i, the gap equation is (Q n — uj n + ifi) 

S(p,w„) _1 = ry • p + i~fiU) n + m + 'Sip, Q n ) , (1) 

I— — oo 

^^l^S{q,U}{)—Y v {q,Lb u p,Lj n ), (2) 

where: ui n = (2n + 1)ttT is the fermion Matsubara fre- 
quency; Q n i = uj n — uji ; is the drcsscd-gluon propa- 
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gator; and T v is the dressed-quark-gluon vertex. (As we 
employ an ultraviolet-finite model, renormalisation is un- 
necessary and m = in Eq. (1) defines the chiral limit.) 
The gap equation's solution can be expressed as 

5(p,w n ) _1 =17 -pA(p 2 ,ui 2 ) 

+i l4 u n C(p 2 ,Cj 2 n ) + B(p 2 ,Cj 2 n ), (3) 

with, e.g., B(p 2 ,w 2 )* = B{p 2 ,u 2 _ n _ 1 ). The dressed- 
gluon propagator has the form 



g 2 D fll ,(k, n n i) = p^D T (k 2 ,n 2 nl ) + p^D L (k'- 



(4) 



where P^ L are, respectively, transverse and longitudinal 
projection operators. Whilst for T ^ ^ /i it is generally 
true that Dt ^ Dli there are indications [26] that for 
T < 0.2 GcV, the domain with which we are concerned, 
it is a good approximation to treat Dt = Dt ='■ Dq. 
For the in-vacuum interaction we use a simplified form 
of that in Rcf. [27]; viz., with n = k 2 



O 2 



, x 47T 2 „ ,2 

D (k) = D—n<T K l a 

V (7 



(5) 



The parameters in Eq. (5) are D and a but they are not 
independent: a change in D can be compensated by an 
alteration of a [24]. For a G [0.3,0.5] GeV, using Eq. (6) 
below, ground-state pseudoscalar and vector-meson ob- 
servables are roughly constant if aD = (0.8 GcV) 3 . We 
usually use a = 0.5 GeV. NB. Eq. (5) is used for illustra- 
tive simplicity, not out of necessity. The status of prop- 
agator and vertex studies can be tracked from Ref. [25]. 

The gap equation is complete once the vertex is spec- 
ified. For the meson spectrum it is now possible to use 
any reasonable Ansatz [28]. Herein we compare results 
obtained using the rainbow-truncation: 



I\,(<f,wz,p,w n ) = 7„ , 



(6) 



the leading term in a symmetry-preserving scheme [23], 
with those produced by the Ball-Chiu Ansatz [29, 30]: 



+{pn + «)/ 



(7) 



with (F = A, B, C) 



Z F (q 2 ,uj?,p 2 ,u> 2 n ) = -[Ftf 2 ,u>?)+F(f 2 ,u; 2 n )], (8) 

A /-a -2 x F(q 2 ,u; 2 )-F(p 2 ,u; 2 n ) 

A F {q ,ui,p u n ) = L 2 _ - 2 , (9) 

Ql Pn 

where, defining u = (0,0,0,1), 7J = 7^ - u^ a u a , 
l)l = u ^laU a ; and p n = (p,uj n + in), q t = (q,u t + ifi). 
The comparison is natural because vertices of the type 
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FIG. 1. Panel (a) - Temperature dependence of the chiral- 
symmetry order parameter in Eq. (11). Chiral susceptibility 
computed in the Wigner phase - Panel (b), and in the Nambu 
phase - Panel (c). In all panels: the Ball-Chiu vertex was 
used, Eq. (7); and fj, = 0, D = 0.5 GeV 2 , m = 0. 



in Eq. (6) are widely used in studies of hadron observ- 
ables [24, 25], and the Ball-Chiu Ansatz provides a semi- 
quantitatively accurate representation of lattice-QCD re- 
sults for important terms in at T = = \i [31]. 

We have solved the gap equation formulated above and 
in Fig. 1 depict the T-dependence of a chiral susceptibility 
[18] and a chiral-symmctry order parameter [20] 



X(0,w ) = ^B(6,lo 2 ), 



(10) 



-(qq)° = N C T ]T tr D f -j^L S m=0 (p,UJ n ) . (11) 

n=-oo •* ^ W ' 

For T < T E the behaviour of the order parameter is typ- 
ical of models without long-range correlations in the gap 
equation's kernel [18]. Namely, initially slow evolution 
from its T = value: (qq)° = (-0.258 GeV) 3 , which sig- 
nals chiral symmetry realised in the Nambu mode; i.e., 
dynamically broken chiral symmetry, and this followed 
by a mean-field transition to a phase with chiral symme- 
try restored; i.e., realised in the Wigner mode. 

The lower panels in Fig. 1 show the chiral susceptibil- 
ity of the Wigner and Nambu phases, which correspond 
to gap equation solutions that are, respectively, within 
the domain of attraction of the B = or B ^ solution 
[32] . A phase is unstable in response to fluctuations if the 
susceptibility is negative, but stable and realisable other- 
wise. With [i — 0, one sees the Nambu phase completely 
replaced by the Wigner phase at T = 124 McV. 

This should be contrasted with the behaviour in Fig. 2. 
At T = the order parameter remains constant with in- 
creasing /1 until fjL a = 0.30 GeV, which is the upper bound 
on the domain of analyticity for our gap equation's kernel 
[33]. On a small domain beyond this; viz., [i 6 (/i Q ,/^), 
with pj£ = 0.314 GeV, the order parameter diminishes 
smoothly, an effect that may be denominated a partial 
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FIG. 2. Analogue of Fig. 1, displaying evolution with chemical 
potential at T — 0. 



restoration of chiral symmetry. For fj, > ^ the order pa- 
rameter vanishes so that chiral symmetry is completely 
restored via a first-order transition. 

The lower panels of Fig. 2 provide additional informa- 
tion. For fx < fj™, with fj™ = 0.286 GeV, the Wigner 
phase is unstable. However, that changes at ^ , when 
the chiral susceptibility in the Wigner phase switches sign 
and thereafter, on the domain < \i < fj,^ , both 
the Wigner- and Nambu-phasc susceptibilities are pos- 
itive. This is the domain of phase coexistence, with 



a metastable Wigner phase for /x^ < fi < [i a and a 
metastable Nambu phase for fi a < /i < [if? . The pressure 
of the phases is equal at fi = fi a ; and the Nambu phase 
is completely displaced by the Wigner phase for fi > fi^ . 
Notably, with an Ansatz for the dressed-quark gluon ver- 
tex, the diagrammatic content of the gap equation's ker- 
nel is generally unknown. However, owing to the insights 
provided in Ref. [34], one can draw these conclusions de- 
spite being unable to calculate an explicit expression for 
the thermodynamic pressure. 

At the onset of the coexistence domain we expect pock- 
ets of deconfined, chirally symmetric quark matter to 
appear in the confining Nambu medium. Their num- 
ber and average volume will increase with fi. The oppo- 
site situation occurs at the termination of the domain; 
i.e., it is the Nambu phase which exists only in pockets. 
For [i e (/i a , /i^ ) , which is the domain of Nambu-phasc 
metastability, the properties of observed hadrons will be 
affected by the partial restoration of chiral symmetry. 

We performed computations at many (/i, T)-values and 
therefrom drew the phase diagrams in Fig. 3. The upper 
panel was obtained with the rainbow truncation, Eq. (6), 
in which case the diagrammatic content of the gap equa- 
tion's kernel is known and one can thus compute the 
dressed-quark component of the pressure. It is given 
by the auxiliary-field effective action evaluated at its ex- 
trema [17]. Within the domain (fi < ^ E ,T > T E ) the 
system exhibits a mean-field transition, which is signalled 
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FIG. 3. Chiral-limit phase diagram in the temperature/che- 
mical-potential plane for strongly-interacting quarks. The 
critical endpoint (CEP) is marked explicitly. Upper panel - 
rainbow vertex, Eq. (6), with D = 1.0 GeV 2 . Lower panel - 
Ball-Chiu vertex, Eq. (7), with D = 0.5 GeV 2 . 



both by: equality of the Wigner- and Nambu-phase pres- 
sures; and coincident singularities in the Wigner- and 
Nambu-phase chiral susceptibilities. The curve track- 
ing the singularity location in the Wigner and Nambu 
susceptibilities bifurcates at the CEP, with a domain of 
phase coexistence opening. Naturally, the curve of equal 
thermodynamic pressure lies within this domain. 

The lower panel in Fig. 3 was obtained using the 
dressed vertex, Eq. (7). Its features are similar to those 
displayed in the upper panel but the domain of phase co- 
existence is smaller. Here, as mentioned above, one can- 
not derive an expression for the dressed-quark pressure. 
Thus it is only with an appreciation of the information 
contained in the chiral susceptibilities that it is possible 
to draw a phase diagram at all. 

It has long been conjectured that confinement is ex- 
pressed in the analytic structure of the dressed-quark 
propagator [17, 24, 25]. Measured in this way, it is signif- 
icant that the models we have considered are members 
of a class in which chiral symmetry restoration is ac- 
companied by a coincident dressed-quark deconfinement 
transition in the chiral limit. 

In Table I we illustrate the response of the CEP's lo- 
cation to changing the vertex or the parameters. Defin- 
ing a confinement length-scale r a = 1/er, it is apparent 
that the CEP rotates toward the temperature axis as 
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TABLE I. Parameter- and vertex-dependence of the critical 
cndpoint and coexistence region. Ac is the width of the co- 
existence region on the /i-axis. Each parameter set produces 
similar results for the in-vacuum values of the so-called vac- 
uum quark condensate (qq)° and the pion's leptonic decay 
constant f n [28]. Lattice-QCD suggests (/x B ,T B )/T c = (1.0- 
1.4,0.93). (Dimensioned quantities in GeV.) 



model 


result 


vertex 


D l/2 


a 


Tc 


Ac 


(^E,r E )/r c 




BC 


0.71 


0.50 


0.124 


0.026 


(1.13,0.89) 


1.27 


BC 


0.71 


0.45 


0.142 


0.075 


(0.85, 0.88) 


0.96 


Bare 


1.00 


0.50 


0.133 


0.220 


(0.98, 0.90) 


1.08 


Bare 


1.02 


0.45 


0.138 


0.280 


(0.80,0.88) 


0.90 



r a is increased. The extreme case is r a = oo, which was 
computed in Rcf. [21] and reported above: (/z E , T E )/T C = 
(0, 1). Models of the NJL-type, as they have been used in 
the current context, represent the opposite limiting case: 
they are expressed via a gap equation in which the con- 
finement length-scale vanishes. From this perspective, it 
is unsurprising that they produce a CEP whose angular 
separation from the /i-axis is significantly smaller. 

We described a method, based on the chiral suscepti- 
bility, which enables one to draw a phase diagram in the 



chcmical-potential/tcmpcraturc plane for quarks whose 
interactions are described by any sensibly-constructed 
gap equation. Thus, in attempting to chart the phase 
structure of QCD using the methods of continuum quan- 
tum field theory, one is no longer restricted to the sim- 
plest class of mean-field kernels: sophisticated quark- 
gluon vertices can be used. The method is general and 
potentially useful in all branches of physics that explore 
the properties of dense fermionic systems. 

A class of models that successfully describes in-vacuum 
properties of it- and p-mesons, exhibits a critical endpoint 
(CEP) in the neighbourhood (/i E , T E ) - (1.0,0.9)T C . 
The CEP's angular separation from the temperature axis 
is a measure of the confinement length-scale: the separa- 
tion decreases as the confinement length-scale increases. 
Furthermore, a domain of phase coexistence opens at the 
CEP. It's size depends on the structure of the gap equa- 
tion's kernel but, other aspects being equal, it increases 
in area as the confinement length-scale increases. Our 
results suggest that illumination of the CEP and its con- 
sequences is within the reach of modern colliders. 
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